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ABSTRACT 

The cosmic far-infrared background (CFIRB) is expected to be generated by faint, dusty star- 
forming galaxies during the peak epoch of galaxy formation. The anisotropy power spectrum 
of the CFIRB captures the spatial distribution of these galaxies in dark matter halos and the 
spatial distribution of dark matter halos in the large-scale structure. Existing halo models of 
CFIRB anisotropy power spectrum are either incomplete or lead to halo model parameters that 
are inconsistent with the galaxy distribution selected at other wavelengths. Here wc present 
a conditional luminosity function approach to describe the far-IR bright galaxies. We model 
the 250 firn luminosity function and its evolution with redshift and model-fit the CFIRB power 
spectrum at 250 fim measured by the Herschel Space Observatory. We introduce a redshift 
dependent duty-cycle parameter so that we are able to estimate the typical duration of the 
dusty star formation process in the dark matter halos as a function of redshifts. We find the 
duty cycle of galaxies contributing to the far-IR background is 0.3 to 0.5 with a dusty star- 
formation phase lasting for ~ 0.3 — 1.6 Gyrs. This result confirms the general expectation that 
the far-IR background is dominated by star-forming galaxies in an extended phases, not bright 
starbursts that are driven by galaxy mergers and last ~ 10 — 100 Myrs. The halo occupation 
number for satellite galaxies has a power-law slope that is close to unity over < z < 4. Wc find 
that the minimum halo mass for dusty, star-forming galaxies with L250 > 10 10 L Q is 2 x 10 11 
M Q and 3 x 10 10 M© at z = 1 and 2, respectively. Integrating over the galaxy population 
with L250 > 10 9 L Q , we find that the cosmic density of dust residing in the dusty, star- forming 
galaxies responsible for the background anisotropies fWst ~ 3 x 10 -6 to 2 x 10 -5 , relative to 
the critical density of the Universe. 



1. INTRODUCTION 

The total intensity of the cosmic far-infrared back- 
ground (CFI RB) is now establ is hed with absolute 
photometry (iPuget et al.l 119961: iFixen et all 119981: 
IDwek et al.l H998L This background originates from 
the UV and optical emission of young stars, absorbed 
by the dust in galaxies and then re-emitted in the 
infrared (IR) wavelengths. Deep surveys with in- 
struments aboard the Herschel Space Observatory 
have started to resolve this background intensity be- 
tween 100 and 5 00 urn to discrete galaxies based on 
resolved counts (1 liver et al.l 120101 : IClements et al.l 
120101 iBerta et al.l 1201 lh . Unfortunately even the 
deepest images of the far-IR sky using PACS and 
SPIRE are limited by source confusion. For exam- 
ple, at 250, 350, and 500 fim, only 15, 10 and 6% of 
the total background intensity is resolved to individ- 
ual galaxies, respectively. Instead of individual de- 
tections, the fainter galaxies responsible for the bulk 
of the CFIRB intensity is studied with st atistics such 
as P (D), the probability of deflections (|Glenn et al.l 
I2010D . and P(k), the angular power spectrum of 
CFIRB anisotro pies resulting from the correlated 
confusion noise (IHaiman fc Knox I [2000l: E nox et aLl 
20011 : IScott fc White I 119991 : INegrello et al.l I2007t 
Amb ard fc Cooravll2007t ). 

While attempts were made to detect the power 
spectrum of CFIRB with Spitzer-MIPS at 160 ^m 
(Lag ache et al.1 l2007t ). and a limited low signal- 
to-noise ratio detection with BLAST (De vlin et al.1 



2009]) in iViero et al.l <|2009f) . the first clear detec- 
tion of the CFIRB anisotropy power spectrum from 
30 arcseconds to 30 arcminute angular scales came 
from Herschel-SPIKE at 250, 350 and 500 (im 
(jAmblard et al.l 120111) . This was soon followed by 
Planck measurements of the CFIRB power spec- 
trum from 5 arcminute to degree angular scales 
(jPlanck Collaboration] 1201 lh . At the Ion ger mm- 
wavelengths, clustering of dusty galaxies can also be 
studied as part of the CMB secondary anisotropy 
studies, where a combination of signals contribute 
to the total power spect rum (|Addison et al.l 120111 : 
lArchidiacono et all [2012). While the arcminute- 
scale ground-based CMB experiments and Planck can 
study the large-angular correlations in the CFIRB 
at linea r scales, the angul ar resolution of Herschel- 
SPIKE (jGriffin et al.ll2010D is such that the measure- 
ments probe the non-linear scales and capture im- 
portant information on how the dusty star-forming 
galaxies are distributed in the dark matter halos. 

First predictions on the CFIRB anisotropy power 
spectrum concentrated on t he linear power spec- 
trum scaled by a bias f actor (IScott fe White"! Il999: 
IHaiman fe Knox I 120001: iKnox et al.l 120011) . Since 
those early studies a popular approach to describe the 
large-scale structure galaxy distribution is to connect 
galaxies to the un derlying dark matter ha lo distribu- 
tion (see review in iCoorav fc Sheth1l2002T) . This halo 
modeling allows a way to describe the galaxy cluster- 
ing power spectrum and correlation function through 
the halo occupation number describing the number 



of galaxies in a given dark matter halo as a func- 
tion of the halo mass. Recent improvements in the 
halo model involve an occupation number description 
that takes into account for the luminosity dependence 
of the satellites through what are now called condi- 
tional luminosity functions (CLFs; lYang et al.1 [2004; 
ICoorav fc Milosavlievicll2005t lCooravll2006al) . 

While an attempt was made to incorporate 
CLFs to describe the CFIRB power spectrum 
(lAmbard fc Cooravl 120071: see recent works i n 
IShang et al.ll201lHWang etaI1l20Tlt IXia et al.ll2012h . 
this was based on a phenomenological model for the 
number co unts and luminosity f unctions (LFs) of far- 
IR sources (jLagache ct al. 2003). The number counts 
and LF measurements from the Herschel Space Ob- 
servatory now allow us to both improve the model 
and extract parameters of the underlying CLF de- 
scription. 

Separately, modeling of recent measurements of 
the CFIRB anisotropy power spectrum with Herschel 
and Planck, and the dusty galaxy signal in CMB 
secondary anisotropy data, is somewhat controver- 
sial. The best-fit pa rameters of the original study 
(jAmblard et al.l 120 111) either had a power-law slope 
for satellites that was steeper than 1.3 or had a re- 
lation between the satellite mass scale M sat and the 
minimum halo mass to host a galaxy M m { n such that 
M,at ~ (3 — 4)M m i n . The galaxy clustering mea- 
surements in the optical band sh ow that the power- 
law slope is slightly l ess than 1 (jZehavi et al.l [2004; 



lAbazaiian et all 120051). while M S a t ~ (1 5 - 20)M m 
(IGao et all 12004 Kravtsov et al.l 120031; iZheng et al.1 
I200a Hansen et all 120091; IShang et al.H201lD . Thels^ 
sue is not limited to the Herschel power spectrum 
since similar conclusions can also be reached with 
fits to the Planck CFIRB power spectrum. Prior to 
Herschel, the low signal-t o-noise CFIRB po wer spec- 
trum reported by BLAST (jViero et al.ll2009D required 
a halo profile that extends out to ~ 4r v i r to fit the 
data, leading to an overestimate of the mean density 
of dark matter in the universe relative to the value in 
f] m that normalizes the dark matter halo mass func- 
tion. 

While a power-law description of the CFIRB power 
spectrum out to £ ~ 2000 is likely adequate for 
the dusty galaxy power spectrum in ground-based 
arcm inute-scale CMB anisotropy data ( Addison ct al. 
120111) . a clear departure from the power-law was 
detected in the Herschel measurement out to I > 
10 4 , indicating the transition between the 2-halo 
and 1-halo term of galaxy clustering. A proper de- 
scription of the Herschel CFIRB power spectrum 
must then move beyond the p ower-law fit to th e 
data. We also re fe r the reader to iPenin et al.1 (|2012[) ; 
iBethermin et~aLl ([20121) for more recent modeling 
of CFIRB, concentrating on the Planc k -meas ured 
CFIRB power spectrum, and IXia et al.l (|2012D for 
modeling of both Planck and Herschel. 

This paper is organized as follows: in the next Sec- 
tion we outline the Herschel data used for this analy- 
sis. In Section 3 we present a revised CLF model for 
the CFIRB anisotropy power spectrum. In Section 4 
we present our results and conclude with a summary 
in Section 5. Throughout this paper we assume the 



fiducial cosmology for the AC DM model of WMAP-7 
results (jKomatsu et al.1120111 ). 

2. DATA USED FOR THE ANALYSIS 

The CFIRB angular power spectrum used f or thi s 
analysis is the same as that of lAmblard et al.l (|2011[ ). 
taken f rom the Herschel Multi-tiered Extra-galactic 
survey (jOliver et al.ll2010l) wi th the Spectral an d Pho- 
tometric Imaging Receiver (jGriffm et al.l I2010D on- 
board the Herschel Space Observatorv (jPilbratt et al"] 
I2010T ) . While the measurements were reported for 
three wavelengths, we concentrate on the 250 /jm an- 
gular power spectrum since it has the highest signal- 
to-noise and the best resolution. 

The luminosity f unction data measur ed by Her- 
schel are taken from lVaccari et aTl (|2010l ) at low red- 
shifts (z < 0.2) at 250/im. The high-z luminosity 
function data extend i ng up to z = 4 at 2 50 /ur n 
are from lEales et all (120101) and ILapi et al.l (|201lD . 
We use data from lEales et al.l (|2010D in 2 bins at 
0.2 < z < 0.4 and 0.4 < z < 0.8. These luminos- 
ity functions are based out of optical and near-IR 
photometric or spectroscopic redshifts for 250 /j,m- 
detected galaxies in the GOODS-North field. To ex- 
tend the 250 //m LF s to higher r edshif ts we make use 
of the results from ILapi et al.1 (|2011f h These lumi- 
nosity functions are somewhat uncertain as they are 
based on the sub- mm photometric redshifts, which 
for each galax y could have an e rror of at least 0.3 
in Az/(1 + z) ([Harris et al.ll2012| ). In any case some 
of that uncertainty is captured by the errors of the 
LF. In future with more exact LFs our model can be 
further improved. 

The low-z luminosity function data are shown in 
Fig.[T]left, while the CFIRB angular power spectrum 
data are in Fig. [2] left. Fig. [3] left shows the high 
redshift luminosity functions data. 

3. HALO-MODEL AND LUMINOSITY FUNCTION 
FORMALISM 

The CLF model we use here to analyze the LF 
and P(k) data is largely ba sed on the model of 
iGiavalisco fc Dickinsonl (|2001l ) and iLee et all (j2009f ). 
One of the main advantages of this improved descrip- 
tion of a halo model is that it clearly connects the 
luminosity function to clustering of galaxies, allow- 
ing to simultaneously constrain the model parame- 
ters using both these observables. The connection be- 
tween 1-point (LF) and 2-point (w(6), P(k)) is based 
on an explicit model of the galaxy luminosity-halo 
mass relation as a function of redshift. When com- 
pared to the standard halo model with galaxy statis- 
tics described by an occupation number, such a lu- 
minosity based approach is capable of accounting for 
the fact that the luminous galaxies are more likely 
to be in more massive h alos. Moreover, the CLF- 
based model description of lLee et alj (|2009f ) is general 
enough to reproduce a wide range of shapes for the 
galaxy luminosity-halo mass relation and its scatter. 
This is advantageous as the shape of this relation is 
expected to be different at far-IR wavelengths when 
compared to the same data at optical wavelengths. 

The fundamental ingredients in this revised CLF 
model are the mass functions for halos and sub-halos 
and the galaxy luminosity-halo mass relation and its 



evolution. The probability density for a halo or a sub- 
halo of mass M to host a galaxy with luminosity L is 
modeled as a normal distribution with 



P(L\M) = 



f/DC 



1no L {M) 



exp 



2a L (M) 2 



,(1) 



where t^dc is the duty cycle factor related to the 
duration of the star formation in the halos (and is 
< ?7dc < !)• More precisely, the duty cycle repre- 
sents a measure of the duration of the star formation, 
£sf, relative to the time interval probed by the survey 
or the observations, At. As discussed in ILee et ahl 
(2009), the ratio tsF/At determines the number of 
halos that can host a detectable galaxy and is hence 
related to the ratio between the number densities 
of galaxies and available halos to host such galaxies 
n g /rih- This is precisely the duty cycle 7/dc we have 
introduced above. 

Note that this description of the duty cycle is dif- 
ferent from the "duty cycle" reported bv lShang et al.l 
(|2011l) for dusty, starforming galaxies in their 
halo/CLF modeling of the Planck CFIRB power spec- 
tra. In their work, the duty cycle is derived by com- 
paring the measured shot-noise amplitude, that dom- 
inates anisotropy power spectrum at small angular 
scales, to a prediction of the expected shot-noise given 
the number density of halos and the observed counts. 
Given that the shot-noise is J dSS 2 dn/dS, where S 
is the flux density and dn/dS is the number counts, 
the shot-noise quoted in their paper is weighted more 
toward s the bright, r are so urces. The model compar- 
ison bv lShang et al.l (|201 ID suggests a duty cycle that 
is close to one suggesting that the CFIRB anisitropies 
are dominated by normal quiescent galaxies. Here we 
provide a precise estimate of the duty cycle down to 
a specific luminosity and as a function of redshift. 

In eq.[TJ a{M) is the scatter in the luminosity-mass 
relation. In this description the scatter can be re- 
lated to the nature of the starformation. High values 
for the scatter with respect to the mean luminosity 
L(M) imply a star formation dominated by starbursts 
while low values for scatter, suggesting a fixed rela- 
tion between halo mass and luminosity, are typical of 
q uiescent, steady star formation (see also discussion 
in ILee et al.ll2009ft . 

The relation between the halo mass and the aver- 
age luminosity L{M) is expected to be an increasing 
function of the mass with a characte r istic m ass scale 
M i and we can write (see ILee et al.l (j2009f )) 



L(M) = L (—J exp 



M 



-0i 



, (2) 



and the scatter can be parameterized in a similar way 



a(M) = cr 



/ M 



\M 0l 



exp 



/ M 



\M 0s 



(3) 



As already discussed bv lLee et al.1 (|2009l ) these pa- 
rameterizations don't have a specific physical motiva- 
tion (except for the requirement of being increasing 
function of mass) , but offer the advantage to explore 
a large range of possible shapes. We need to con- 
sider the total halo mass function, that is the number 



density of halos or sub-halos of mass M. The con- 
tribution of halos nh(M) is taken to be t he Sheth & 
Tormen relation (jSheth fc Tormenl 119991 ) . The sub- 
halos term can be modeled through the number of 
sub-halos of mass m inside a parent halo of mass M p , 
N(m\M p ). The total mass function is then 

n T (M) = n h (M) + n sh (M), (4) 

where n s h(M) is the sub- halo mass function 



(5) 



n sh (M)= N(M\M p )n h (M p )dM p . 



We pa rameterize N{m\M) as in Ivan de Bosch et al.l 
(2005) 



N(m\M) 
where 

7 = 



7 / m 



pT{l-a sh ) \Mp sh 

fsh 



exp 



r(i - a sh , i/p sh ) - r(i - a sh , io- 4 //U) 



Mp sh 



(7) 



Here T(x) is the incomplete gamma function and 
fsh is the sub-hal o mas s fraction. As shown in 
Ivan de Bosch et aTl ((2005), where the model is cali- 
brated using numerical simulations, both the normal- 
ization and the slope of the sub-halo mass function are 
not universal and depend on the ratio between the 
parent halo mass and the non linear mass scale, M*, 
defined as the mass scale where the rms of the den- 
sity field cr(M, z) is equal to the critical over-density 
required for the spherical collapse S c (z) . The term 
f s h in eq. [7] is fitted by the relation 

\og[(f sh )} = [0.4(log(M/M„) + 5)] 1/2 + 2.74 , (8) 

in numerical simulations and we make use of it in this 
study. 

The best-fit rel ation for the slope param eters a s h 
and fish found by Ivan de Bosch et al.l (12005) is 



a sh = 0.966 - 0.028 log(M/M*) , 



(9) 



and fi s h — 0.13, independent of M. With this de- 
scription the total number of free parameters in the 
CLF model is 9 with 4 parameters for the luminosity- 
mass relation in eq. [2] 4 parameters for the scatter in 
eq. [3l and the duty cycle parameter ?7dc • 

If the same luminosity-mass relation applies 
to both halos and sub-halos, then the product 
P(L\M)riT(M)dLdM gives the number densities of 
galaxies with luminosity L in halos or sub-halos of 
mass M. The luminosity function is then 



cf>(L)dL = dL / dMP{L\M)n T {M) . 



(10) 



The formalism introduced above also allows us to 
construct the halo occupation distribution (HOD) in 
a simple way. The contribution of central galaxies is 
simply the integration of P(L\M) over all luminosities 
above a certain threshold Lq either fixed by the survey 
or a priori selected so that 



(N C (M)) 



L>L,, 



P(L\M)dL, (11) 



-t 



which, in absence of scatter, reduces to a step function 
9(M — Mo) as expected. For the satellite galaxies, 
the HOD is related to the sub-halos 



spectrum is 



P 2,l (k)- 



(N S (M)) 



L>L„ 



dL / dmN(m\M)P(L\r^2) 



dM(N T (M))u(k, M)n h (M)b(M$ 8) 



xPlin(fc) , 



(19) 



The total HOD is then 



where Pn n (fc) is the linear power spectrum and 
b(M) is the lin e ar bi as factor calculated as in 
(N tot (M)) L>Lmin = (N h (M)) L>Lmin + (N sh (M)) L>Lr il3 p 00 ^ fc Shet S (EQQ2). The total galaxy power 

' spectrum is then P g (k) = P lh (k) + P 2h (k). 

As the observations are anisotropics on the sky 
projected along the line of sight, the observed an- 
gular power spectrum can be related to the three- 
dimensional galaxy power spectrum thr ough a red- 
shift integration along the line of sight (|Knox et al.l 
l200l : 



The model described so far holds at a given redshift. 
The duty cycle parameter tjdc and the luminosity- 
mass relation with its scatter are expected to have a 
redshift evolution. Here we are attempting to fit LF 
data at a variety of redshift bins between < z < 4. 
To account for the redshift evolution of the parame- 
ters we assume that the parameters that the describe 
the low redshift (z < 0.2) dusty galaxy population 
are different from those for the high redshift galaxies. 
Moreover, the high-z data extend from z = 0.2 to 
z = 4 and are divided in to 6 redshift bins. We thus 
fit a total of 7 duty cycle parameters, one for each 
of the bins and do not attempt to constrain the duty 
cycle variation with a parameterized approach on the 
redshift evolution. For the evolution of the galaxy 
250 /im luminosity-halo mass relation we account for 
the possible redshift evolution by introducing another 
parameter, pm and rewriting the mass scale Mo/ as 



M Ql (z) = M 0lz<0 . 2 {l + zf 



(14) 



where we allow the evolution to follow the assumed 
power-law form. 

Once the HOD is defined, it is possible to cal- 
culate the one-halo and two-halo terms of the far- 
IR anisotropy power spectrum. First, we define the 
power spectrum in terms of redshift-dependent three- 
dimensional clustering and will later project them 
along the line of sight to calculate the angular power 
spectrum of CFIRB anisotropies. Here we assume 
that the central galaxy is at the center of the halo 
and that the halo radial profile of satellite galaxies 
within dark matter halos follow that of the dark mat- 
ter give n by the Navarro, F renk and White (NFW) 
profile (|N avarro eFaIIH99l . The one-halo term is 
then 



1 



P lh {k) = — / dM{N T (N T - l))u{k,MYn h We) 
ni 



where u(k, M) is the NFW profile in Fourier space 
and n g is the galaxy number density 



n g = / dM(N g (M))n h (M). 



(16) 



The second moment of the HOD that appear in eq. [15] 
can be simplified as 

(N T (N T - 1)) ~ (N T ) 2 - (N h ) 2 , (17) 

and the power index p for the NFW profile is p = 1 
when (Nt(Nt — 1)) < 1 and p = 2 otherwise 
(|Lee et al.ll2009D . The two- halo term of galaxy power 



(~) j,A:)j,A :l/'„(//.\ . -^20) 



where x is the comoving radial distance, a is the scale 
factor and j v (z) is the mean emissivity at the fre- 
quency v and redshift z per comoving unit volume 
that can be obtained from the LFs as 



j v (z) = I dL(j)(L,z)— . 



(21) 



This model does not rely on the assumption of a num- 
ber counts shape or an evolution. We are able to di- 
rectly model-fit the mean emissivity as a function of 
redshift. 



ai 0.22 ±0.10 

ft 0.70 ±0.05 
log(M oi /M ) 11.5 ±1.7 

log(L /i Q ) 9.6 ±2.4 

Vdc (0 < z < 0.2) 0.54 ±0.26 

TABLE 1 

Best-fit parameter values and their 1 — a uncertainties 

FROM THE LOW- REDSHIFT LUMINOSITY F UNCTION DATA OF 
IVACCARI ET AL.I (|2010h 
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0.57 ±0.03 






0.19 ±0.02 


\og{L' /L Q ) 




9.58 ±0.02 


C*250 




0.78 ±0.16 


IDC (0.2 < 2 


< 0.4) 


0.43 ±0.07 


IDC (0.4 < z 


< 0.8) 


0.30 ±0.04 


■que (1.2 < z 


< 1.6) 


0.16 ±0.01 


r)DC (1-6 < 2 < 2.0) 


0.19 ±0.01 


r^Dc (2.0 < 2 


< 2.4) 


0.33 ±0.01 


t)DC (2.4 < 2 


< 4.0) 


0.31 ±0.02 


PM 




-4.32 ±0.09 



TABLE 2 

Best-fit parameter values and their 1 — a uncertainties 

FROM THE COMBINATION OF ANGULAR CFIRB POWER 
SPECTRUM AT 250/iM AND HIGH-Z LUMINOSITY FUNCTION DATA. 



4. RESULTS AND DISCUSSION 

In the revised CLF model outlined above, in princi- 
ple, we have 23 free parameters: 7 duty cycle param- 
eters, 8 parameters for the luminosity-mass relation 
and its scatter at low redshifts, 7 parameters for the 
same relations at high redshifts plus the parameter 
Pm for the (1 + z) redshift evolution of the mass scale 
(eq.|H). 

Separately, the CFIRB power spectrum contains 
the contribution from Galactic cirrus, in addition 
to the extragalactic anisotropies tra c ed by the faint, 
dusty galaxies. In lAmblard et al.1 (|2011[ ) the au- 
thors accounted for this contamination assuming the 
same cirrus power-law power -spectrum from measure- 
ments of IRAS and MIPS (jLagache et all 12007ft at 
100/im and extending it to h igher wavelengths us - 
ing the spectral dependence of lSchlegel et al.l ( |1998|) . 
Such a frequency scaling resulted in an overesti- 
mated cirrus correction, as noted by the Planck 
team (jPlanck Collaboration] 1201 ID in their analy- 
sis of the CFIRB power spectrum compared to the 
Herschel power spectrum. This is primarily due to 
the fact that the c irrus is likely overestimated in 
ISchlegel et afl (fl998h as IRAS 100 (im also contains 
the extragalactic background intensity. To avoid bi- 
asing our power spectrum low by an overestimated 
cirrus cor rection, we re-fi t the raw power spectrum 
data from lAmblard et al.1 (|201lD . Here we adopt the 
same power-law cirrus fluct uation power spectrum 
used in lAmblard et al.1 (|2011l ). with P( k) oc fc " with 
n = —2.89 ± 0.22 as measured by (|Lagache et al.1 
2007]) . However, we rescale the amplitude of the cir- 
rus power spectrum with a dimensionless factor C250 
that we keep as a free parameter and model-fit that 
as part of the global halo model. This implies another 
free parameter in our model, leading to a total of 24 
parameters. Given the large volume of the parameter 
space, a MCMC analysis (see below) through the full 
parameter space is very time-consuming. Moreover, 
it is unlikely that the information carried by the cur- 
rent data is able to constrain such a large number of 
free parameters. 

We hence simplify the analysis as follows. We first 
fit the low-redshift parameters to the z < 0.2 250 11m 
lumin osity function measurements from lVaccari et all 
( 2010) by only varying the 4 parameters related to the 
luminosity-halo mass relation and the duty cycle at 
low redshift, r]nc(z < 0.2). We assume no scatter in 
the luminosity-mass relation. When fitting to high 
redshift data the total number of free parameters is 
11: 3 parameters for the L — M relation (aj, /3 Z ', L' ol , 
while Mo2 iZ <o.2 is kept fixed to the value found for the 
z < 0.2 LF), plus the power index pu to account for 
the evolution of M i , 6 duty cycle parameters and 1 
amplitude for the cirrus contamination, 6*250. When 
model fitting to the measured angular power spec- 
trum at 250 11m we calculate the total theoretical Ct 
as the sum of the Ci for the low-z HOD found in the 
previous fit and the CV calculated at z > 0.2. This 
allows us to treat two types of galaxy populations 
that are contributing to the Herschel galaxy popu- 
latio, the low-z (z < 0.1) dust in late- type galaxies 
and the dusty spher oidal galaxies at high redshifts 
(|Lagache et al.l |2003|) . and to account for a possible 



redshift evolution of the others L-M relation param- 
eter (L , ai, /?;). Here we consider galaxies brighter 
than L > 5 x 10 7 L Q for the low redshift model, while 
we use L > 10 9 L Q to model-fit the high redshift data. 
These values are consistent with the flux cut of the 
galaxy samples considered. In order to account for 
the uncertainty in the exact value of L m in we have ver- 
ified that an order of magnitude change in the value 
of L m in leads to changes in the power spectra of the 
order of 5 to 6% which is comparable to the lrr er- 
ror bars of the data. We find that the exact value of 
Lmin, within an order of magnitude, does not change 
the results considerably. 

To model-fit the data we implement a Markov 
Chain Monte Carlo (MCMC ) analysis using a mod- 
ified version of the cosmoMC iLewis fe Bridle I (2002) 
package. The results for the low-z z < 0.2 250 (im 
luminosity function data are shown in Fig. [1] and in 
Table [TJ where we show and tabulate the best-fit to 
the luminosity function data and the best-fit values 
for the CLF parameters involved with the LF descrip- 
tion, respectively. In Fig. Q] we show the HOD calcu- 
lated for the best-fit values of the parameters and its 
uncertainties. 

The best-fit model to the angular power spectrum 
data and high-z luminosity functions is shown in Fig- 
ures [5JDH In Table [5] we tabulate the best-fit param- 
eter values (where the prime is used to distinguish 
the parameters for the high-z luminosity-mass rela- 
tion from those for the low-z case) and in Fig. [4] we 
show the probability contours for the luminosity-mass 
relation parameters. 

These results show that the model is able to fit 
the data even assuming no scatter in the luminosity 
mass relation. The HOD shows a sharp cut-off at a 
mass of about log(M min /M ) ~10.8atz = 0. This 
quantity could be compared to the threshold mass 
of the standard hal o model and it is in agreement 
with the results of lAmblard et al.l (|201lD . where it 
was found that log(Af m in/-M©) — 11.5 with a simple 
HOD for the dusty galaxies. However it should be 
noted that in this work we are not fitting directly 
the value of M m i n as it is not a free parameter in 
our model. Thus a d irect comparison of our work to 
lAmblard et al.1 (|2011l ) may not be appropriate. 

Both the values ment ioned above are d ifferent 
from the recent results of iShang et al.1 ()201lD . where 
the authors used an improved version of the halo 
model including a luminosity-mass relation to ana- 
lyze t he Planck-based CFIRB anis otropy power spec- 
trum (|Planck Collaboration J|2TjTTI ). They found that 
the most efficient halo mass scale for starformation 
is M c ff ~ 10 12 ' 65 Mq, which is closer to the typical 
value of optical galaxie s in the standard halo model 
(jCoorav fc Sheth H2002tlAbazaiian et al.l|2005ft. Nev- 
erthel ess, again, as already noted in IShang et al.l 
( 201 j]), the model used there is different from the one 
used in lAmblard et al.1 (|2011| ) (and from the one used 
in this paper) , and a direct comparison between M e f 
and M m in may not be accurat e. A proper compa rison 
of our model to the results of IShang et al.1 (|2011l ) can 
be done through the effective halo mass scale: 

f Nt 
M eS = / dMn h (M)M-^- , . (22) 
J ng 
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Fig. 1. — Left: 250 (im luminosity function data at z < 0.2 and the best-fit model. Right: The halo occupation number at z < 0.2. 
The minimum halo mass to host a galaxy with luminosity L > 5 ■ 1O 7 L0 is log(M m i n /Af0) ~ 10.8 and the power-law slope of the 
satellite galaxies with halo mass is ~ 0.98. The grey region represents the 68% confidence level and the three lines at the top right 
of the figure correspond to power-law slopes of 0.5, 1.0 and 1.5. 
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Fig. 2. — Left: Best-fit model description to the angular power spectr um data at 250/xm. W hen fitting to the angular power 
spectrum data we removed the shot- noise contribution as determined by Amblard et al. (2011). Dashed and dotted lines are the 
1-halo and 2-halo term, respectively. The shaded region is the 68% confidence level galactic cirrus contribution, as determined by 
the free parameter C250 included in the model-fit. Right: The halo occupation number for the high-z galaxies. The HOD slope 
at the high mass-end is ~ 0.96. The grey region represents 68% confidence level and the three lines correspond correspond to 
power-law slopes of 0.5, 1.0 and 1.5. 




Fig. 3. — Left: Best-fit model description to the high-redshift luminosity functions from Ealcs ct al. (2010) and Lapi ct al. (2011). 
Right: Satellite fraction at each redshift bin. 

Integrating over our HODs, we find log 10 (M c ff) = 12.63 at z = 0, which is comparable to the re- 
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Fig. 4. — 68% and 95% confidence level contours for the parameters of the luminosity-mass relation L'q-{5\ {left) and L'^-a^ [right). 
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Fig. 5. — Duty cycle ?7dc as a function of redshift. Over the 
redshift range of 1 < z < 4, rjrjc ~ 0-2 to 0.4. 




Fig. 6. — The 250 /im luminosity-halo mass relation and the 
68% confidence level region at z = 0, 1 and 2. 



suits of S hang et al.l f|201 lh ■ It is also comparable 
to results for optical g alaxies ()Coorav fc Sheth | [2002; 
lAbazaiian et al.l [20051 . At z = 1.4, where the emis- 
sivity peaks (see Fig. Qand discussion below), we find 
M eff = 11.45. 

The contribution of satellite galaxies in our best- 
fit model becomes efficient at a mass scale M sat ~ 
17M m ; n and the HOD of satellite galaxies has a 
power-law behavior of oc M s with the power-law slope 



Fig. 7. — 250 /^m emissivity predicted using the CLF 
model of this s t udy c ompared to the model prediction of 
[Bethermin et al. <20TlI) . The shaded region correspond the 
the 68% confidence level from the MCMC model fits to the 
measurements used here. 
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Fig. 8. — The cosmic density of dust Qdust vs redshift. The 
CLF model prediction (gray region) calculated with the best- 
fit luminosity functions of Fig. [3| and compared to other works 
(see text). 

of s ~ 0.98. We show in Fig. [5] that this slope re- 
mains close to unity also for the high-z HOD where 



we find s ~ 0.96. The relation M s: 



17M miD be- 



tween minimum halo mass and the halo mass scale 
at which satellites appears and the slope of the 



satellite galaxies occupation number are consistent 
with expectations from numerical si mulations and re- 
sults obtained for optical galaxies (IGao et ah! 12 004; 



Kravtsov et al.ll2003t IZheng et al.ll2005HHansen et al l 
2009 HShang et al.ll2011l) . 

Our result o n the power-law slope with s ~ 0.98 is 
different from lAmblard et al.l ()2011l ) where a higher 
slope (> 1.6) was found. The difference with s ~ 1 
and s > 1.6 between the two works comes in part 
from the fact that the model we have presented here 
accounts for that the brighter satellites are in more 
massive halos and in part from our rescaling of the cir- 
rus amplitude. As discussed i n iPlanck Collaboration.! 
(|2011j) . lAmblard et al.l (|2011f) overestimated the cir- 
rus contamination and underestimates the clustering 
power spectrum of the cosmic infrared background. 
We have refitted the cirrus amplitude as part of the 
joint fit to the power spectrum measurements using 
both a halo model and a power-law power spectrum 
for cirrus. 

The simple halo occupation number used in 
lAmblard et al.l (| 2 1 If ) is not able to make a dis- 
tinction between the luminosity and mass of satel- 
lite galaxies. Moreover, in the current description, 
the shape of the luminosity-mass relation determines 
the strength of the 1-halo term, while in the stan- 
dard halo-model, the strength of the 1-halo term is 
mainly determined by the slope of the halo occupa- 
tion number of satellite galaxies. As the 1-halo term is 
clearly detected in the Herschel CFIRB power spec- 
trum the previous results were biased by an incor- 
rect model that attempted to model-fit high signal- 
to-noise power spectrum measurements. We are also 
finding an amplitude for the cirrus contamination 
which is smaller than 1 : C250 = 0.78±0.16. A reduced 
cirrus contamination requires a reduced relative am- 
plitude between the 1-halo and 2-halo terms that can 
be achieved with a lower slope for the satellite con- 
tributions. We note also that the fit to the cirrus 
amplitude confirms th at the previous analys is of the 
same clustering data in lAmblard ctldl (|201 If) overes- 
timated the cirrus contribution as a l ready has been 
found by the IPlanck Collaboration] ()201lD , where a 
factor ~ 2 of difference was found in the cirrus con- 
tamination. The model used in this work is hence able 
to alleviate the tension between Planck and Herschel 
analysis. 

In terms of the duty cycle parameters, at z < 0.2, 
where the constrain came only from the luminosity 
function data, we found a weak constraint on the 
duty cycle, that is DC = 0.54 ± 0.48 at 95% con- 
fidence level. The amplitude of the luminosity func- 
tion is in fact affected by both the duty cycle and 
other parameters of the luminosity mass relation, in 
particular Lq and Mqi. Thus, in the absence of other 
constraints, the duty cycle can not be measured ef- 
ficiently. On the other hand, the parameters of the 
L — M relation also determine the HOD and the rel- 
ative amplitude of 1-halo and 2-halo terms. Combin- 
ing clustering measurements and luminosity function 
data can hence strongly improve the constraints on 
the duty cycle. Also the z < 0.2 luminosity function 
is better determined in narrow, deeper surveys com- 
pared to the case of the z < 0.2 LF. All this results 



in a determination of high-z duty cycle parameters 
with relatively small uncertainties. We have found 
that the duty cycle generally decrease with redshift 
until z ~ 1 and slowly increase again for higher z 
(see tableland Fig.[S]). The values are in the range 
?7dc ^0.16-0.4 when .2 < z < 0.4. 

In iShang et al.l (|2011l ) the authors have excluded 
very low duty cycle values (?7dc < 0.05) finding that 
the Planck power spectrum data favor duty cycles 
close to unity. Our results lie in the mi ddle between 
these two extreme cases. However in IShang et al.l 
(|2011[ ) the duty cycle parameter is estimated by com- 
paring the shot noise predicted fo r a fixed ring to 
the re sults of the empirical model of IBethermin et al] 
(2011), rather than fitting to the data. Moreover the 
model used here differs in the description of the lu- 
minosity mass-relation and introduces a possible red- 
shift evolution both for the duty cycle and for the 250 
/xm luminosity-halo mass relation. Interestingly our 
constraints are much more similar to the results of 
ILee et all (|2009D . where the model we are using here, 
with some differences in the L-M relation and its red- 
shift dependence, was also used to analyze the UV lu- 
minosity function and two point correlation function 
data of starf orming galaxies i n the range z = 4 — 6. 
In that work ILee et alj (|2009h used ?7dc as an input 
parameter rather than as a free parameters and found 
that extremely short (770c < 0.1) and extremely long 
(jjdc > 0.7) duty cycles are ruled out at the 90% con- 
fidence level. Our results also suggest a mid range for 
t]dc- The agreement could imply that a large fraction 
of the UV-se l ected starforming galaxy sample studied 
in ILee et al.l (|2009|) could also be responsible for the 
CFIRB anisotropics. This could be directly tested 
via a cross-correlation between the two datasets and 
such studies are expected in the near future given 
the Herschel imaging of some of the wide area legacy 
fields. 

In Section 12] we described the duty cycle in terms 
of the duration of the starformation tsF in the halos 
with respect to the time interval At covered by the 
survey. A long duty cycle 770c ~ 1 implies a starfor- 
mation time scale that is tsF 3> At, while small duty 
cycles with rjuc ~ correspond to the opposite case 
with rjgp <§C At. The central value with ?7dc ~ 0.5 im- 
plies tsF — At. In terms of the physical time, once ac- 
counted for the time interval spanned by each redshift 
bin, the duty cycles listed in Table [5] correspond to a 
starformation phase lasting for tgF — 0.3 — 1.6 Gyr. 
Such a long starformtaion t imescale is c onsist ent with 



what has been sug gested inlLapi et al.l ( 20111) and the 
physical model of iGranato et al.l (|200U I2004D . Such 
a long time scale rule out models where the CFIRB 
is dominated by gas-rich mergers with tsF — 10 — 100 
Myr. 

It is worth noticing that in this analysis we are as- 
suming a duty cycle that is independent of mass. It 
could very well be that the duration of starformation 
depends on the halo mass. Given the large number 
of free parameters in the analysis we are not able to 
parametrize a possible mass or luminosity dependence 
of ?7dC) but we regard this possibility as a future im- 
provement to this model. 

Our analysis suggests that the L-M relation has 



a redshift dependence, with Mqi that decreases for 
higher redshifts (see eq. HH) . Decreasing M i is equiv- 
alent to increase the characteristic luminosity of the 
luminosity function. The fit to the high-z luminos- 
ity function data is shown in Fig. [3J In particular 
we find M i oc (1 + z )-^±om_ Although a direct 
comparison is complicated because of the very differ- 
ent models used, we n ote that this re s ult is similar 
to the evolution seen in iLeFloc'h et all (pOOl . where 
the characteristic luminosity has been found to have 

a redshift dependence oc (1 + z) ' -°- 2 . 
The total luminosity-mass relation calculated as: 



more massive (and hence less numerous) halos with 
respect to central galaxies. This result for h a t is also 



L(M) = L(M) + J N(m\M)L(m)dM, (23) 

is shown in Figure |6] The shaded regions represent 
the lcr uncertainty and it can be seen that the data 
are able to constrain the luminosity-mass relation 
with good precision even for the z < 0.2 case using 
only the low redshift luminosity function data. The 
luminosity-mass relations we are finding show a lin- 
ear behavior. However the luminosity functions are 
steep at the low-faint luminosity end. This result in 
a tension when the observed turnover is attempted 
to be expl ained through the abun dance matching ap- 
proach in Bcthcrmi n et al.l (|2012T ) . We consider this 
a natural consequence of the fact that in this work 
we are attempting to fit simultaneously datasets in 
a large range of redshifts together with anisotropy 
power spectrummeasurements. Moreover we observe 
that there is no clear visible turnover in the data 
that we are fitting without imposing any prior or con- 
straint on the faint-end of the luminosity functions. 
The faint-end description, both in data and in mod- 
els, should be further improved. 

In Figure [7] we show the emissivity correspond- 
ing to the best-fit model, calculated according to 
equation (|2T1) and c ompared to the emissiv ity of the 
parametric model of iBethermin et ahl ([20111 ) (see also 
iPenin et al.ll2012fl . The extended tail at z > 3 is due 
to the constant and the fast increase of the luminos- 
ity function with red shift. We have v e rified that us- 
ing the emissivity of [Bethermin et ahl (|2011l ) implies 
a few percent difference in the best-fit values of the 
CLF parameters in the model presented here, com- 
parable to the lcr error bars. Future analysis may 
require however a different redshift parameterization 
of the average luminosity-mass relation. 

We show also the satellite fr action for the four 
high r edshift bins calculated as Ivan de Bosch et ahl 
(200j): 



fsat(L) 



1 



dMP(L\M)n T {M) , (24) 



M' 



where M' is the mass scale where there is one galaxy 
brighter than L. The satellite fraction is an impor- 
tant test for galaxy formation models and to establish 
the properties of galaxy-halo relation. We find that 
the fraction is about 22 — 25% at L = 1O 9 L0 and 
decreases quickly to less than 5% at 10 11 L Q while we 
don't find a significant redshift dependence. The de- 
creasing behavior with mass is due to the fact that 
satellite galaxies at a given luminosity are located in 



inspect to cenuai f^aitixics. ims icauit 101 j s n t is aisu 

i n agreement with van de Bosch etap (pOM ) : [Cooravl 



(20060): ICounon et all (|20lT^ 



Finally in Fig. [5] we show the fraction of dust with 
respect to the critical density of the Universe p, cal- 
culated as 



^dust — 



1 

P0 



dL(j){L,z)M dust (L), (25) 



where Mdust is the dust mass correspond ing to a 
given IR luminosity and we use Eq. 4 in iFu et all 
(2012|). The results are com pared to those in 
Fig. 7 of iMenard fe Fukugita.l (|2012t) where O dust 
has been determined with reddening of metal-line 
absorbers. In Fig. [8] we also show other esti- 
m ates of the mass dens i ty of dust as summarized 
bv IMenard fc FukueitaTI (l2012fl from iFukueita et al.l 
(200 41): iDriver et al.l (|2007D : IMenard et al.l (|2oTo[ h 



Fuk ugital (1201 If). W e have combined the points from 
Menard ct al. "or the dust contributions of ha- 

los and those from iMenard fc Fukugita. I (|2012[ ) m a 
single set of data points, under the assumption that 
the amount of dust in halos doesn't evolve signifi- 
cantly with redshift. We parameterize the opacity 
with a power low kd oc v^ d with the power index in 
the range (3d = 1.5 — 2. The calculation requires the 
spectral energy distribution of dust and we assume a 
thermal black-body spectrum wit dust temperature 
in the range T — 25 — 35X. We allow for a large 
range in dust temperature, taken as a uniform prior, 
to allow for the rang e of values seen in current data 
(jAmblard et al.ll2010l) . In equation ^ we integrate 
over luminosities L min > 1O 9 L . However in this 
calculation the choice of L m i n is less relevant, since 
the uncertainty on fWst is dominated by the large 
range of temperatures and spectral indices consid- 
ered. The gray region correspond to the prediction 
for these parameters ranges using the best-fit lumi- 
nosity functions of Fig. [3] 

5. CONLCUSION 

We have presented an analysis of the Herschel- 
SPIRE CFIRB power spectrum at 250 pva and the 
luminosity functions up to z — 4. We use a con- 
ditional luminosity function approach to model the 
far-IR bright galaxies. We have modeled the 250 /im 
luminosity function and its evolution with redshift in- 
troducing a redshift dependent duty-cycle parameter. 
This description represents an improved version of the 
halo-model that offers an advantage by accounting for 
the luminosity dependence of the satellite galaxies as 
a function of the halo mass. The underlying ingredi- 
ent is the galaxy luminosity-halo mass relation. 

We have found that current Herschel data are able 
to constrain the model despite the high number of free 
parameters. The results of our analysis indicate that 
the cosmic far-IR background is dominated by star- 
forming galaxies in an extended phase of starforma- 
tion rather than bright starbursts that are fueled by 
gas-rich mergers. We found duty cycles correspond- 
ing to a dusty starformation phase lasting ~ 0.3 — 1.6 
Gyr, which is in agreement with previous analysis of 
starforming UV-selected galaxies at high redshifts. 



1U 



We have also found that the halo occupation num- 
ber for satellite galaxies has a power-law slope that 
is about 0.98 over the redshift range < z < 4. This 
solves the tension between previous analysis of the 
same Herschel power spectrum data and other deter- 
minations of the halo occupation number for galaxies 
in the literature. Finally we have estimated the cos- 
mic density of dust residing in the dusty, starforming 
galaxies responsible for the cosmic far-IR background 



anisotropics to be fidust ^ 3 x 10 6 to2xl0 5 . 
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